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(57) Abstract 

Flow parameters in a fluid delivery assembly are determined by monitoring pressure responses and processing those responses along 
with information regarding the fluid flow. In one aspect, a processor controls the pump to pump flow quantities in accordance with a 
pseudorandom code. Portions of the resulting pressure signal sensed are then decoded in accordance with the pseudorandom code. An 
estimate of the equilibrium pressure is generated from the decoded pressure values, while a summation of the pressure samples is generated 
from the undecoded pressure signals. The resistance to fluid flow of the system is determined from the estimated equilibrium pressure and 
pressure summation. For low flow rates, a processor controls the pump to pump fluid in a series of fluid boluses, with each fluid bolus 
delivered in the beginning of a separate timeslot. The equilibrium pressure is measured at the end of each timeslot, and a summation of the 
pressure samples is generated from the pressure signals. For high flow rates, the pump is controlled to vary the flow rate and the change in 
pressure is divided by the change in flow to directly determine the resistance. A resistance display continuously displays the resistance of 
the system. The pseudorandom coding and decoding can be used to filter out pressure-response crosstalk caused by multiple fluid infiision 
segments feeding into a common line. 
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WO 00/13726 



PCI7US99/18334 
FLUID FLOW RESISTANCE MONITORING SYSTEM 



This is a continuation-in-part of application Serial No. 
08/305,904, filed September 13, 1994. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The invention relates to fluid delivery systems. More 
particularly, the present invention relates to monitoring the 
resistance to fluid flow in a fluid delivery system infusing 
5 fluid to a patient. 

Description of Related Art 

There are a variety of situations where fluid is infused 
to a patient. Applications of fluid delivery systems include 
(but are by no means limited to) intravenous infusion, intra- 

10 arterial infusion, infusion of enteral solutions, infusion of 
medication to the epidural space, and diagnostic infusion to 
determine vascular characteristics of the arterial, urinary, 
lymphatic, or cerebrospinal systems. 

Fluid delivery systems for infusing fluid to a patient 

15 typically include a supply of the fluid to be administered, 
an infusion needle or cannula, an administration set 
connecting the fluid supply to the cannula, and a flow 
control device, such as a positive displacement infusion 
pump. The administration set typically comprises a length of 

20 flexible tubing. The cannula is mounted at the distal end of 
the flexible tubing for insertion into a patient's blood 
vessel or other body location to deliver the fluid infusate 
to the patient. The flow control device often is a 
peristaltic-type pump that acts on the flexible tubing to 

25 force the fluid through the tubing of the administration set 
to the cannula and into the patient. One commonly used flow 
control device is a linear peristaltic type pump having 
several cams and cam-actuated fingers that sequentially 
occlude portions of the flexible tubing along a pumping zone 

30 to create a moving zone of occlusion. 
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During an infusion procedure, events may occur that 
interfere with proper delivery of fluid to the patient, such 
as an occlusion of the administration line. It is desirable 
to detect these events as soon as possible so that they can 
5 be remedied. 

A common technique for detecting such events and for 
evaluating fluid delivery system status is to monitor the 
pressure in the administration set. Variations in pressure 
can indicate problems with fluid delivery. For example, an 
i0 increase in pressure over a selected threshold may indicate 
an occlusion in the system. Similarly, a drop in pressure 
can indicate an empty fluid supply or other fluid delivery 
system fault. 

A problem with determining fluid delivery system status 

15 by monitoring pressure alone is the slow speed at which 

pressure can build when the system is operating at a low flow 
rate. At low flow rates, the energy per unit time introduced 
into the flow path is very small. Accordingly, it may take a 
considerable amount of time for the pressure to build up 

20 enough to exceed a threshold and indicate an occlusion. 
Additionally, with a relatively low pressure threshold, 
patient movements such as coughing, sneezing, and sitting up 
can cause the pressure to momentarily exceed the pressure 
threshold, thus creating a false alarm of a fluid delivery 

25 system fault. Another problem with monitoring pressure alone 
occurs when the delivery cannula becomes mislocated within 
the interstitial tissue matrix, causing a rise in pressure. 
The amount of resulting pressure rise is dependent upon flow 
rate. For example, at a flow rate of 10 ml/hr # the rise in 

30 pressure is typically only about 10 mm Hg; at a flow rate of 
2 ml/hr, the rise in pressure is typically only about 
2 mm Hg. Such small relative changes are difficult to detect 
from instantaneous pressure readouts, or even from pressure 
trends, because of the presence of other sources of change, 

35 such as patient movement, as well as venous pressure, 

elevation of the system components, and the flow rate itself. 
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As has been noted in U.S. Patent 4,898,576 to Philip, 
the measure of the resistive part of the fluid line impedance 
can be used to monitor the condition of the fluid line. One 
technique used in actively monitoring the resistance, rather 
5 than merely waiting for pressure to build up, is the 

alteration of the flow rate. The change in the pressure over 
the change in the flow rate has been found to accurately 
indicate the resistive part of the fluid impedance in the 
system when adequate time is allowed for the pressure to 

10 reach equilibrium at each rate. This technique has been 
found to be effective at higher flow rates with their 
accompanying higher pressures. A change in these higher flow 
rates is accompanied by a rapid and measurable change in 
pressure. Because of the rapid pressure response to the flow 

15 rate changes, the flow rate can be varied about the selected 
flow rate without any significant clinical effect on flow 
uniformity. 

However, at lower flow rates, the clinical requirement 
of flow rate uniformity restricts the magnitude of the 

20 changes to the flow that can be imposed on the fluid line. 
It is thus undesirable to alternate between different flow 
rates to obtain different pressure responses for determining 
resistance due to the detrimental effect on flow uniformity 
the flow changes would have as well as the relatively long 

25 length of time required to obtain those pressure responses. 
Various models of pressure and resistance can allow 
accurate resistance measurements. For example, as described 
in U.S. Patent No. 5,087,245 to Doan, which is incorporated 
herein by reference, a technique for determining flow 

30 resistance which allows for a non-linear relation between 
pressure and flow and a time-varying impedance (resistance 
and compliance) involves inducing a bolus of fluid in the 
infusion system and monitoring the resulting pressure wave 
and the pressure decay response. Injecting a known quantity 

35 of fluid causes a resulting pressure wave that then decays to 
an equilibrium pressure. Using the equilibrium pressure and 
the pressure decay response, fluid resistance can then be 
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determined even when a non- linear relation between flow and 
pressure exists and when the impedance (resistance and 
compliance) are time -varying via the following equation: 

A [(Pit) - P 0 ) dt 

Resistance = = J 

Af f-F(t) dt 



where: / F(t) dt = Q = the known delivered quantity 
5 of fluid, 

P(t) = the change in pressure over time, 

P 0 = the equilibrium or offset pressure, 

Ap = the area under a pressure response 
waveform, and 

10 A f = the area under a fluid flow waveform. 

Thus, by injecting a known quantity of fluid through the 
infusion system, monitoring the resulting pressure as it 
decays to an equilibrium pressure, and determining an 
integral of the difference between the equilibrium pressure 

15 and the pressure response, the resistance to fluid flow can 
be determined. However, after the known quantity of fluid is 
injected, further quantities of fluid can not be injected 
through the system (i.e., further flow steps are not 
initiated) in order for the pressure to decay to the 

20 equilibrium pressure. In some situations, such as where the 
fluid resistance is relatively high, a relatively long pause 
in fluid injection may be necessary to allow the pressure to 
reach equilibrium pressure. Depending on the particular 
application, such long delays between fluid flow steps may be 

25 undesirable. 

As set forth in pending U.S. patent application Serial 
No. 06/305,904, pseudo- random binary sequence (PRBS) codes 
have been used to effectively eliminate the delays in 
reaching equilibrium pressure by creating "virtual" 

30 waveforms. However, due to the high processing requirements 
of PRBS coding and decoding procedures, PRBS codes have only 
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been used in combination with linear and non- time -variant 
models of pressure and resistance, such as the following 
equation: 

[F(t) 

P{t) = Resistance * F( t) + ^ 

Compliance 



5 Where F = Flow rate, 

P = Pressure , and 

Resistance and Compliance are stationary values 
(i.e., values that do not vary with time or flow) . 

Such a linear and non-time variant estimation technique is 

10 relatively accurate over fluid resistances between 0 to 1500 
fluid ohms (where 1 fluid ohm = 1 mm Hg per liter per hour) . 
However, because fluid impedance (i.e., resistance and 
compliance) is in actuality time-variant and non-linear, the 
above-cited estimation technique has reduced accuracy where 

15 the fluid resistance exceeds 1500 fluid ohms. 

Note that fluid resistance is a part of the total system 
fluid impedance. Fluid impedance is a function of the system 
compliance, inertance, and resistance. 

The causes of the non-linear, time-variant resistance 

20 relationship include the viscoelasticity of the flexible 
tubing, which slowly contracts following application of a 
positive pressure transient. Additionally, the biochemical 
and rheological aspects of the patient's fluid system, such 
as the blood flow in the human body, further complicate the 

25 pressure/flow relationship. 

There are several applications of parenteral infusion 
systems where both low flow rates and high resistances are 
encountered. For example, infusing parenteral liquids into 
small children and infants, and particularly into premature 

30 infants, can involve low flow rates and high resistances. 

An additional problem with monitoring fluid infusion 
systems is caused by various "noise" sources that can degrade 
the pressure monitoring. Such noise sources include movement 
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of the fluid infusion system and patient movement such as 
breathing and ambulation. Additionally, using more than one 
pump and administration set to inject fluid through a single 
fluid delivery system cannula can introduce noise that can 
5 interfere with individually monitoring the flow in the 
various administration sets. 

Hence, those skilled in the art have recognized a need 
for a fluid delivery monitoring system that can detect a 
fluid delivery fault condition faster and with improved 

10 specificity than prior systems at low flow rates. There is 
also a recognized need for a system that can compensate for 
the existence of offset pressure while maintaining clinically 
acceptable flow patterns, and that can detect partial or 
"soft" occlusions that may result in pressure changes that 

15 are too small to be noticed through conventional pressure 
monitoring systems. Additionally, it has also been 
recognized that there is a need for a system that is less 
sensitive to other sources of pressure changes in the conduit 
such as those caused by other pumps on the same fluid line. 

20 It is further desirable to have a system that is accurate 
over a wide range of resistances and that is less sensitive 
to noise effects. The present invention satisfies these 
needs and others, 

SUMMARY OF THE INVENTION 

25 Briefly and in general terms, the present invention is 

directed to a system that monitors one or more flow 
parameters in a fluid delivery assembly. In a fluid delivery 
system in which a flow control device acts on a fluid conduit 
to control the movement of fluid through the conduit, the 

30 system for monitoring one or more flow parameters comprises a 
pressure sensor coupled to the conduit for providing pressure 
signals in response to the pressure sensed in the conduit, 
and a processor that receives the pressure signals, processes 
those pressure signals, and determines a flow parameter based 

3 5 on said processing. 
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In one more detailed aspect, the processor determines 
flow resistance. 

In a further aspect, the processor controls the flow 
control device to cause various flow rates to exist in the 
5 conduit, with the application of specific flow rates 
dependent upon whether the selected flow rate is high, 
medium, or low. The processor applies different resistance 
measurement techniques depending on the selected flow rate. 
In a further aspect, a pseudorandom code is used to 

10 create a repeating, non-uniform flow pattern, which induces 
pressure responses which are measured and then decoded to 
compute the estimated equilibrium pressure or total offset 
pressure. The pressure responses over a pseudorandom code 
period are also summed. The pressure response summation and 

15 the estimated equilibrium pressure are then used to determine 
the flow resistance. 

In yet another aspect, where the selected flow rate 
falls bel ow a low threshold, the processor determines a 
timeslot length and controls the flow control device to 

20 deliver a bolus of fluid at the beginning of the timeslot. 
The pressure response is monitored to determine an 
equilibrium pressure, and the sum of the pressure response is 
determined. The resistance is calculated using the 
determined equilibrium pressure and the pressure response 

2 5 sum . 

In a further more detailed aspect, where the selected 
flow rate exceeds a high threshold, the processor controls 
the flow control device to cause a plurality of different 
flow rates to exist in the conduit. The processor then 

30 processes the difference in the pressures and the difference 
in the flow rates to determine the impedance to flow. 

In yet a further aspect, the processor controls the flow 
control device to deliver a closely spaced series, or 
"trill", of flow waveforms at the beginning of a timeslot. 

35 In another aspect that may be used with a fluid delivery 

system that includes two or more fluid infusion segments, 
each of which may include a separate fluid source and a 
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separate flow control device acting on a separate fluid line; 
each of the two or more fluid infusion segments feeds into 
common fluid line that delivers fluid to a patient; and where 
at least one of the fluid infusion segments includes a 
5 processor that controls the flow control device, the 

processor uses a pseudorandom coding and decoding process to 
filter pressure-response crosstalk caused by the other fluid 
infusion segments. 

In another more detailed aspect of the invention, the 
10 system determines an estimate of signal quality and noise. 

Other features and advantages of the present invention 
will become more apparent from the following detailed 
description of the invention when taken in conjunction with 
the accompanying drawings . 

'15 BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a simplified block diagram of a system for 
detecting abnormalities in a fluid line incorporating the 
principles of the invention as applied to an intravascular 
fluid infusion system. 
20 FIG. 2 is a graphical representation of fluid flow and 

corresponding pressure responses . 

FIG. 3 is a graphical representation l of different fluid 
flow waveforms and corresponding pressure responses. 

FIG. 4 is a simplified block diagram illustrating the 
25 process for generating an estimated pressure waveform using 
PRBS coding and decoding. 

FIG. 5 is a simplified flowchart showing a process for 
determining resistance according to a preferred embodiment of 
the invention. 

30 6 is a graphical representation depicting sample 

motor step and Superstep assignments. 

FIG. 7 is a graphical representation of PRBS code 
generation. 

FIG. 8 is a graphical representation of pressure 
35 responses showing baseline sample averaging regions. 
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FIG. 9 is a graphical representation of PRBS decoding. 
FIG. 10 is a graphical representation of decoded 
pressure responses showing block averages. 

FIG. 11 is a simplified block diagram showing median 
5 filtering. 

FIG. 12 is a graphical representation of flow waveforms 
resulting from adjacent pump steps. 

FIG. 13a is a graphical representation of flow waveforms 
according to an acceleration/deceleration waveform. 
10 FIG. 13b is a graphical representation of flow waveforms 

according to trilled acceleration/deceleration waveforms. 

FIG. 14 is a graphical representation of a pressure 
waveform and corresponding timeslot sum array. 

FIG. 15 is a simplified flowchart showing a process for 
15 determining resistance according to one embodiment of the 
invention . 

FIG. 16 is a graphical representation of pressure 
waveforms at very low flow rates. 

FIG. 17 is a simplified diagram showing a system for 
20 monitoring resistance in a fluid infusion system having two 
fluid infusion segments. 



DETAILED DESCRIP TION OF THE PREFERRED EMBODIMENTS 
Referring now to the drawings with more particularity, 
wherein like reference numerals in the separate views 

2 5 indicate like or corresponding elements, there is shown in 
FIG. 1 a block diagram of a resistance monitoring system 10 
incorporating aspects of the current invention. The 
resistance monitoring system is coupled to the conduit 12 of 
a fluid delivery system in which a pump motor 14 drives a 

30 pumping mechanism 16, which in the embodiment shown comprises 
a rotating cam shaft 18 coupled to the pump motor 14 and 
moving a series of peristaltic elements 20. The peristaltic 
elements 20 operate on the conduit 12 to move fluid from a 
fluid source 22, through the conduit 12, and into a patient 

35 24 via a cannula 26. 
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A user input device 28, such as a keypad, provides 
operator instructions, such as flow rate selection, to a 
processor 30. The processor 30 controls the operation of the 
pump motor 14 driving the pumping mechanism 16. A motor 
5 position sensor 32 determines the position of the motor 14 
and pumping mechanism 16 and provides a position signal to 
the processor 30. 

Located downstream from the pumping mechanism is a 
pressure sensor 34 coupled to the conduit 12 to sense 

10 pressure in the conduit. An analog-to-digital converter 36 
( w A-to-D") receives the analog pressure output signals from 
the sensor 34 and converts them to a digital format at a 
particular sample rate controlled by the processor 30. The 
processor 30 receives the digital pressure signals, processes 

15 them as described in more detail below and calculates the 
resistance to flow. A display 38 presents the resistance. 
One or more alarms 4 0 are provided to indicate an 
unsatisfactory resistance level. As a backup to the 
resistance monitoring system, the alarm system 4 0 may also be 

2 0 activated when the pressure exceeds a set threshold. Where 
such a backup system is present, the processor 30 compares 
monitored pressure values against a pressure threshold. If 
the monitored pressure (or, to prevent false alarms, an 
average of several recent monitored pressure values) exceeds 

25 the pressure threshold, the alarm system 40 will be 

activated. Such a pressure threshold will typically be 
relatively high to prevent false alarms. In a preferred 
embodiment, the pressure threshold is 600 mm Hg. 

The selection of a flow rate is made at the keypad 28 

30 and is received by the processor 30. The user may also 

select at the keypad 28 a high or a low resistance range. In 
one embodiment, the high range is from 0 to 2000 fluid ohms, 
and the high range is from 0 to 6000 fluid ohms. The 
selection of range is used to drive the display and alarm 

35 systems, as well as to determine resistance percent (as set 
forth below regarding equation 4) . The user selects the high 
or low resistance depending on the particular application. 
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In a preferred embodiment, the processor 30, after 
receiving the selected flow rate from the user input 28, 
compares 4 2 the selected flow rate against flow rate 
thresholds to determine if the selection is a "high" flow 
5 rate, "medium" flow rate, or a "low" flow rate. Two 

thresholds, one low and one high, are used to divide the 
high, medium, and low flow rates. A selected flow rate 
exceeding the high threshold is considered high, a selected 
flow rate falling below the low threshold is considered low, 

10 and a selected flow rate between the thresholds is considered 
a medium threshold. In a preferred embodiment, the low 
threshold is 0.5 milliliters per hour, and the high flow rate 
is 50 milliliters per hour. 

As one aspect of the resistance monitoring system shown 

15 in FIG. 1, three different approaches to determining the 

fluid system resistance are used so that a wide range of flow 
rates and resistances may be provided by the fluid delivery 
system with continuous, accurate resistance determination. 

For high flow rates, a bi-rate approach is used, whereby 

20 the processor 30 selects 44 two or more different flow rates 
at which the motor will run. The flow rate commands are 
provided to the motor controller 46 which in turn causes the 
pump motor 14 to act on the conduit 12 through the pump 
mechanism 16 to pump fluid through the conduit 12 at those 

25 discrete rates. The rates are selected in one embodiment to 
average to the selected flow rate for the purpose of 
maintaining flow uniformity. 

The pressure responses to the flow waveforms in the 
conduit 12 are monitored by the pressure sensor 34 with 

30 digital pressure signals provided to the processor 30 by the 
A-to-D converter 36. Those pressure signals are processed 
48, which for the high rate involves full pump revolution 
averages -for the last revolution in each bi-rate range. The 
output is then used with flow table information 28 to 

35 calculate 50 the fluid resistance. The calculated resistance 
may be filtered 52, such as through a moving average or 
moving median filter, to improve the accuracy of the 
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resistance calculation. The filtered resistance is shown in 
the display 38 and sent to the alarm system 40. 

Varying the flow rate at relatively high flow rates 
causes a large and rapid pressure response as discussed 
5 above, and the resistance to fluid flow of the system can be 
determined relatively rapidly in accordance with the 
following: 



R = 



?2 ~ Px 
^2 " ^1 



where: R = resistance 
10 F x = first flow rate 

F 2 = second flow rate 
P x =* pressure at the first flow rate 
P 2 = pressure at the second flow rate 

Another relationship that may be used in determining a 
15 resistance percentage and that takes into account the sample 
rate is: 

j°0 _ k»0 



J?% = (SCALE) 



M N 



*M - Flo 



where: R% = resistance in percent; 

SCALE = a scale factor, equal to (1/20), or 
20 (100/(2000 fluid ohms)), for a low user- 

selected resistance scale of 0 to 2000 
fluid ohms, and equal to (1/60), or 
(100/(6000 fluid ohms)), for a high user 
selected resistance scale of 0 to 6000 
25 fluid ohms; 

P = pressure in mm Hg; 

F = flow in liters per hour; 
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M = samples in one revolution at the high rate; and 
N = samples in one revolution at the low rate . 

Therefore, for high flow rates (which in the embodiment 
shown are above 50 ml/hr) , a u bi-rate" approach is used where 
5 the flow rate is varied at two or more flow rates about the 
selected flow rate. The responsive pressure signals are 
monitored. In this approach, the changes in pressure 
resulting from the flow rates are used to directly calculate 
resistance. The flow rates selected and the length of time 

10 that they are each applied are based on averaging to the 

selected flow rate, so that there is no significant clinical 
effect by altering the flow rates. 

Further details regarding determination of fluid 
resistance at high flow rates using the above-cited high- low 

15 flow rate technique are discussed in detail in pending U.S. 
patent application Serial No. 08/305,904, which is 
incorporated herein by reference. 

For low flow rates, which in one embodiment are flow 
rates below 0.5 ml/hour, resistance determination is 

20 performed using a second technique. The drive controller 

does not use PRBS encoding. Instead, the user selected flow 
rate is used by the motor controller 46 to divide the pump 
cycle into a number of timeslots TS, with all timeslots of 
equal length, as shown in FIG. 2a. Timeslot length is 

25 maximized, and the motor is controlled to deliver a fluid 
volume Q TS in each timeslot, preferably as a single bolus 65 
toward the beginning of the timeslot. The resulting pressure 
waveform 66 is shown in FIG. 2b. By maximizing timeslot 
length, sufficient time is permitted for the pressure in the 

30 conduit to decay to the equilibrium pressure P c , which, to 

reduce the impact of noise and other spurious signals, can be 
calculated as the average of the last portion 67 (i.e., the 
"tail") of the waveform in a particular timeslot. 

Referring again to FIG. 1, the pressure is monitored by 

35 the sensor 34, the pressure signals are A/D converted 36, and 
the tail portion of the pressure signals from a timeslot is 
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used to calculate a signal average 54, such as a Timeslot 
Tail Average (TTA) . The TTA value is used to determine P 0 
56. The A/D converted pressure signals are also sent to 
signal processing to generate a summation of the pressure 
5 signal values. The pressure summation and P D value are used 
to calculate the resistance 50. The resistance value may be 
filtered 58 to increase accuracy and remove noise effects. 
The filtered resistance value is then provided to the display 
38 and alarm system 40. 

10 For medium flow rates, which in one embodiment are flow 

rates between 0.5 ml/hr and 50 ml/hr, the actual flow rate is 
varied about the selected flow rate in accordance with a 
pseudorandom pattern of variation about the selected rate. 
The pseudorandom pattern is provided by a pseudorandom 

15 encoder, which in the embodiment shown is a PseudoRandom 

Binary Sequence (PRBS) encoder 60, the operation of which is 
discussed in greater detail with respect to FIG. 7. The PRBS 
code is provided to the motor control 46, which uses volume 
data (from lookup tables 62), the PRBS code, and the user- 

20 selected flow rate to determine a desired repeating, non- 
uniform flow pattern, as discussed in greater detail below 
with respect to FIGS. 6 and 7. The motor control 46 also 
divides the pump cycle into a number of timeslots TS. The 
motor control 46 then causes the pump motor 14 to act on the 

25 conduit 12 through the pump mechanism 16 to pump fluid 
through the conduit 12 at the determined repeating, non- 
uniform flow pattern. 

FIGS, 2c and 2d show sample flow and pressure waveforms 
from a PRBS code. In the example of FIG. 2c, a PRBS code of 

3 0 length 7 and having elements oflllOlOOis used to 
control fluid infusion. One element (i.e., a 1 or a 0) of 
the PRBS code is assigned to each timeslot. For each 
timeslot having a corresponding PRBS element of 1, a single 
small fluid bolus of volume Q TC is injected, as shown in FIG. 

35 2c, preferably at the beginning of the timeslot. The 
resulting pressure waveform 68 is shown in FIG. 2d. 
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Referring again to FIG . 1, for medium flow rates, the 
resulting pressure in the conduit 12 is monitored by a 
pressure sensor 34 that outputs pressure signals. All or 
some of the pressure signals are averaged 54 to create 
5 average pressure values, such as BaseLine Sample Averages 
(BLSAs) . A Baseline Sample Average is an average of a 
portion of the pressure signals in a timeslot, such as the 
tail portions 69 depicted in FIG. 2d. 

Referring again to FIG . 1, the BLSAs are decoded 64 in 

10 accordance with the pseudorandom pattern of variation, and 
the decoded average pressure values are used for estimated 
equilibrium pressure (P D ) determination 56. 

Note that pressure signal averaging 54 is included to 
decrease the processing requirements for decoding and to 

15 eliminate noise/spurious signals. The system could also 
operate by decoding the pressure signals without any 
averaging. However, decoding is processor- intensive , as set 
forth below with respect to FIG. 9. Accordingly, decoding 
all the pressure signals, instead of just average pressure 

20 values, could substantially increase the processing 
requirements . 

The calculated P D values may be filtered 49, such as 
through the median filter process described below with 
respect to FIG. 11, in order to reduce noise effects and 

2 5 increase the accuracy of the P c value. The filtered P c value 
is provided for resistance calculation 50. 

The A/D converted pressure signals are processed 48 to 
provide a pressure summation value. The pressure summation 
value and filtered equilibrium pressure P c are then used to 

30 calculate the resistance 50. The calculated resistance may 
be filtered 58, such as through the median filter process 
described below with respect to FIG. 11, in order to increase 
the accuracy of the resistance value. The filtered 
resistance is provided to the display 38 and alarm system 40. 

35 Note that the P 0 and resistance values can be updated 

with each new timeslot, and does not require an entirely new 
PRBS code cycle to be completed. For example, as shown in 
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FIG. 2d, updated P G and resistance values can be determined 
over the range 70a, which encompasses the first M timeslots 
shown (TS a through TS g , where codelength M = 7) . Upon 
receiving data from another timeslot, such as TS h , updated P Q 
5 and resistance values can be determined over the new range 
. 70b, which encompasses timeslots TS b through TS h . Upon 
receiving data from still another timeslot, such as TS Lt 
updated P Q and resistance values can be determined over the 
new range 70c, which encompasses timeslots TS C through TSj. 

10 Thus, even though M timeslots are used to determine P 0 and 

resistance, the system does not have to wait M timeslots for 
an updated calculation of P Q and resistance. Instead, each 
time new pressure date is received for a new timeslot, an 
updated calculation of P Q and resistance calculation is 

15 performed using data from the new timeslot and the most 
recent (M - 1) timeslots. 

The rationale for using the pseudorandom pattern is 
outlined in FIG. 3, which shows various flow waveforms and 
their corresponding pressure waveforms. In FIG. 3a, an 

20 induced flow waveform 71a comprises a single bolus of fluid 
72a injected at the beginning of a time period T. The bolus 
72a has a volume Q T . The resulting pressure waveform 74a has 
an initial peak 76a that decays to an equilibrium pressure P D 
at or near the end of the time period T. The equilibrium 

25 pressure can be used to accurately determine the flow 

resistance in the fluid delivery system. However, waiting 
for the pressure to decay to the equilibrium pressure can 
take a long period of time, such as a minute or more. 

In many circumstances, the settling time required for 

30 the pressure to decay to the equilibrium pressure is too long 
for practical use in intravenous fluid delivery systems. 
Additionally, delivering a large, single bolus of fluid can 
be undesirable. It is often preferable to deliver a series 
of smaller boluses 72b of fluid equally spaced over the time 

3 5 period T, as depicted in the flow waveform 71b in FIG. 3b. 
However, such equally spaced small boluses 72b can result in 
a pressure waveform 74b such as that shown in FIG. 3b, where 
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there is insufficient time between flow boluses 72b for the 
pressure to decay to the equilibrium pressure. In such a 
case, the equilibrium pressure cannot be directly measured. 
In one aspect of the current invention, a series of 
5 small boluses of fluid are delivered over the time period T, 
with the small boluses spaced over a number M of timeslots 
using a pseudorandom code whose length M equals the number of 
timeslots into which the period T is divided. In the example 
shown in FIG. 3c, The pseudorandom sequence is a PseudoRandom 
10 Binary Sequence (PRBS) whose length M equals 7. The PRBS 
code shown is 1110100, where 1 indicates a bolus of 
fluid delivered in the particular timeslot TS, and 0 
indicates that no bolus of fluid is delivered in that 
timeslot TS. 

15 As shown in FIG. 3c, the pressure waveform 74c resulting 

from the PRBS-coded fluid delivery shows that there is not 
sufficient time between boluses 72c for the pressure waveform 
74c to decay to an equilibrium pressure. Accordingly, the 
equilibrium pressure cannot be directly measured. However, 

20 by decoding the entire pressure waveform of FIG. 3c using the 
same PRBS code (1110100) in a decoding process, an estimated 
pressure waveform 74d can be generated, as shown in FIG. 3d. 
The estimated pressure waveform 74d corresponds to a virtual 
single bolus of fluid 72d delivered at the beginning of the 

25 time period T. Moreover, the final portion of the estimated 
pressure waveform 74d provides a close approximation P Q est of 
the equilibrium pressure P Q that would have actually occurred 
if sufficient time were allowed between flow boluses for the 
pressure to properly decay. 

30 Accordingly, even though the actual delivered flow and 

the actual measured pressure follow the waveforms 71c, 74c 
set forth in FIG. 3c, by decoding the measured pressure 
waveform* 74c with the PRBS code, a "virtual" pressure 
waveform 74d such as that set forth in FIG. 3d can be 

35 calculated, with the "virtual" pressure waveform 74d 

corresponding to the pressure waveform 74a from FIG. 3a that 
is desired for P D determination purposes. 



WO 00/13726 PCT/US99/18334 

-18- 

However, decoding the entire pressure waveform, or even 
substantial portions thereof, requires decoding hundreds of 
individual pressure signals. Because decoding is a 
processor- intensive procedure, it is desirable to reduce the 
5 number of signals decoded. It has been found that by taking 
an average pressure value 78 from each timeslot TS, as shown 
in FIG. 3e, and decoding just one average pressure value for 
each timeslot, the resulting estimated pressure waveform 74f , 
shown in FIG. 3f, can be used to determine equilibrium 

10 pressure, and the resulting estimated equilibrium pressure 
P Q est is a close approximation of the estimated equilibrium 
pressure that would have been generated by decoding all of 
the pressure signals without averaging. 

Accordingly, even though the actual delivered flow and 

15 the actual measured pressure follow the waveforms 71c, 74c 
set forth in FIG. 3c, by averaging the pressure signals and 
then decoding the averaged pressure values with the PRBS 
code, a "virtual" pressure waveform 74 f such as that set 
forth in FIG. 3f can be calculated, with the "virtual" 

20 pressure waveform corresponding to the pressure waveforms 
74a, 74d from FIGS. 3a and 3d that are desired to determine 
the equilibrium pressure. 

FIG. 4 further illustrates the process of creating a 
pseudorandom fluid flow, monitoring the resulting pressure, 

25 decoding the resulting pressure to create a "virtual" 

pressure, and determining an estimated equilibrium pressure 
P 0 . At 80, a desired flow rate of Q/T, where Q = a selected 
quantity of fluid and T = the time period, is shown. Rather 
than deliver the amount Q as a single bolus of fluid in the 

30 time period T, a PRBS coder derives a PRBS code 82, which in 
the embodiment shown is 1 l l o 1 0 0. The motor control 
divides the bolus of fluid into four separate but 
approximately equal boluses of volume Q TS which are delivered 
in accordance with the pseudorandom pattern of the PRBS code 

35 over the time period, as shown at 84. A transducer or other 
pressure sensor 86 produces pressure signals to measure the 
resulting pressure 88. 
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An averaging process 90, such as block averaging, is 
used to average a portion of the pressure signals from each 
timeslot, as depicted at 92. In the embodiment depicted, 
only the tail portion 94 of each timeslot is averaged. The 
5 resulting averaged pressure values are fed into a PRBS 
decoder 96. 

The PRBS decoder 96 processes the averaged pressure 
values and determines "virtual" pressure values 98. These 
"virtual" pressure values are used to determine equilibrium 
10 pressure (P 0 ) , which is then used to calculate the fluid 
resistance . 

FIG. 5 illustrates in part the process for determining 
fluid resistance using the pseudorandom flow pattern coder 
and decoder to determine equilibrium pressure. The user 

15 selects a desired flow rate 28 for delivery to the patient 24 
through the conduit 12. A PRBS encoder/modulator 60 
determines a desired PRBS code and timeslot length based on 
the user-selected flow rate, and the motor control 46 causes 
the pump motor 14 and mechanism 16 to induce flow in the 

20 conduit 12 in accordance with the PRBS code and timeslot 
length . 

The sensor 34 monitors the resulting pressures in the 
conduit and provides pressure signals, which pass through the 
A/D convertor 36 for determination of the fluid resistance. 

2 5 The pressure signals are used to determine a pressure 

summation 48 over a number of timeslots equal to the PRBS 
codelength. Portions of the pressure signals are also 
averaged 54 to determine tail averages or Baseline Sample 
Averages (BLSAs) for individual timeslots. The BLSAs are 

30 then demodulated 64, using an inverse of the PRBS code. The 
demodulated BLSAs are used for equilibrium pressure P Q 
determination 56. 

In calculating the equilibrium P c values for various 
timeslots, the estimated values of P D can vary in quality and 

35 accuracy due to assorted factors, including system noise, 

large offset changes, and hard occlusions. These factors can 
create large changes in baseline and/or impedance to occur, 
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thus causing the decoded output to become unstable, which in 
turn causes the estimated P D value to have decreased quality 
and accuracy. To prevent gross miscalculation of resistance 
due to such conditions, the system determines signal quality 
5 102 by computing a Normalized Sum Absolute Delta's (NSAD) 
parameter, which provides an estimate of noise. If the NSAD 
exceeds a set NSAD threshold for a selected timeslot, the 
estimated P Q calculated for the timeslot is discarded and 
replaced with the P Q from the most recent timeslot to have a 

10 NSAD below the NSAD threshold. Further details of NSAD 
calculation and use are discussed in greater detail in 
connection with FIG. 11. 

In addition to calculating the NSAD, the accuracy and 
quality of the estimated P Q values may be further reinforced 

15 by filtering each P Q estimate through a P D median filter 104. 
Details of the median filter are discussed in greater detail 
below with respect to FIG. 11. 

The estimated P 0 value and pressure summations are used 
for fluid resistance calculation 50 via an equation such as 

20 the following: 

Resistance - ^ - -*■> dt 

A 



f (Fit) dt 



where: / f (t) dt = Q = the known delivered quantity 

of fluid over the PRBS codelength, and 

P 0 = the estimated equilibrium or offset 
pressure . 

25 The estimated resistance can be updated for each 

timeslot, with the calculation using the pressure summation 
and estimated P 0 value from the current timeslot plus the 
pressure summations of the previous M-l timeslots, so that 
data from M timeslots are used to estimate the resistance. 

3 0 This use of data from multiple timeslots increases the 

accuracy and quality of the updated resistance estimate. The 
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accuracy and quality of the updated resistance can be further 
improved by using a resistance median filter 58, similar to 
the median filter 104 described for use with the equilibrium 
pressure P c . The median filter for the resistance stores the 
5 most recent resistance estimates in an array and determines a 
filtered resistance value. Details of median filter 
operation are discussed below in greater detail with respect 
to FIG. 11. In the embodiment shown in FIG. 5, the filtered 
resistance value is output to the warning/alarm system 4 0 as 

10 well as to a display 38, such as a visual bar graph display 
and/or a written trend plotter. 

Referring now to FIGS. 6 through 13, specific elements 
of the medium flow rate system and method are described in 
greater detail, beginning with flow delivery. Typical fluid 

15 infusion pumps are driven by stepped motors having a number 
of steps per pump revolution. Due to the nature of the 
pumps, the volume delivered by different steps can vary 
widely. Some pump steps may even deliver negative volumes 
(i.e., negative flow). To provide consistent delivered fluid 

20 volumes, in a preferred embodiment of the invention the pump 
cycle is divided into a plurality of Supersteps, with each 
Superstep delivering an approximately equal volume of fluid 
Qss. 

The combination of motor steps into supersteps may be 
25 dynamically determined using updated information regarding 
time -variant or pressure- variant changes in fluid flow per 
individual steps. However, in a preferred embodiment, 
division of steps into supersteps is performed using look-up 
tables specific to a particular pump type. 
30 In the embodiment shown in FIG. 6, the pump motor 

provides 200 steps per revolution, with the steps being 
indexed from 0 to 199. Over one pump revolution, the pump 
delivers a volume of fluid Qrev, typically expressed in 
microliters. However, the volume per step (i.e., Qs) varies 

m 

35 widely from step to step, with some steps even providing 

negative volumes. To equalize the flow, the individual motor 
steps are grouped into a smaller number of movements called 
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Supersteps, with the volume of each Superstep (i.e., Qss) 
substantially the same as the volume of other Supersteps. In 
the embodiment shown in FIG. 6, there are 32 Supersteps per 
pump revolution. 

5 Because the volume pumped per motor step varies widely 

between steps, the number of steps in individual Supersteps 
must vary to provide Supersteps with approximately equal 
fluid volumes. For example, in the grouping shown in FIG. 6 
(which is shown only as an example), Superstep 0 (i.e., Qss 0) 

10 includes just two motor steps (Qs 0 and Qs^ , Superstep 1 
(i.e., QssJ includes three motor steps (Qs 2 , Qs 3 and Qs 4 ) , 
and Superstep 29 (i.e., Qss 29 ) includes 101 motor steps (Qs 88 
to Qs 18e ) . Even though different Supersteps consist of widely 
varying numbers of motor steps, each Superstep delivers 

15 approximately the same fluid volume Qss as other Supersteps. 

In a preferred embodiment of the invention, one or more 
superstep volumes (Qss) are delivered in each timeslot in 
which fluid is delivered (e.g., each timeslot represented by 
a binary code of 1) . The number of Supersteps of fluid 

20 delivered in each fluid-delivery timeslot varies according to 
the selected flow rate and other system parameters. (No 
fluid volumes are delivered in non- fluid-delivery timeslots, 
e.g., timeslots with a PRBS code bit of 0.) Representative 
numbers of Supersteps per "flow" (i.e., fluid-delivery) 

25 timeslot are shown below in Table A. 

Timeslot length is determined as a function of the 
codelength M, the user selected rate, and the number of 
codes/revolution at the selected rate. Timeslot length is 
preferably computed to ensure that the specified number of 

30 codes are executed per revolution and that the volume 

produced per code divided by the total time per code produces 
the user- selected flow rate. In a preferred embodiment, the 
codelength M, and codes/revolution values are specified by a 
table lookup as a function of the selected rate, and the 

35 timeslot length is either calculated or specified as a 

function of the selected rate. In a preferred embodiment, 
timeslots are selected to result in the period T being 
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approximately one minute, thus ensuring that the average flow 
rate over a period of no more than T will precisely meet the 
requirements of the operator. 

Table A sets forth values of codelength M and 
5 codes/revolution in one embodiment of the invention. In the 
table values set forth in Table A r the values were selected 
to ensure that the duration of the code is no less than 
approximately 52 seconds over the flow rate range 0.6 to 49.9 
ml/hr for a revolution volume of 183 microliters. Note that 

10 the values set forth in Table A are specific to a particular 
pump apparatus, notably the IVAC Signature Edition, Pump 
(versions 7100 and 7200) , for which the invention divided the 
200 pump steps into 32 Supersteps. However, other pump 
motors and mechanisms are also applicable to the invention, 

15 and other values of M, BIC, etc. may be selected for various 
pumps and other parameters. 



TABLE A 



Selected 
Flow Rate 
(ml/hr) 


Codelength 
M 


Codes Per 
Revolution 


Supersteps 
Per "Flow" 
Timeslot 


0.1 - 0.5 


n/a 


n/a 


1 


0.6 - 1.4 


7 


8 


1 


1.5 - 3.0 


15 


4 


1 


3.1 - 6.1 


15 


2 


2 


6.2 - 12.4 


15 


1 


4 


12.5 - 24.9 


15 


1/2 


8 


25.0 - 49.9 


31 


1/4 


8 



The PRBS codes are preferably of a length that assures 
that they will evenly fit into a complete pump cycle. To 

30 assure an even fit in the pump cycle, the PRBS code should 
have a number N of flow-delivery timeslots (i.e., timeslots 
during which a superstep of fluid is delivered, or "active" 
timeslots) , such that the total number of Supersteps per pump 
cycle (which in the example cited is 32) can be evenly 

35 divided by N. In the example shown in Table A, a PRBS code 
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of length 7 will have 4 flow-delivery timeslots, during each 
of which one Superstep of fluid will be pumped into the 
conduit by the pump. Since there are 4 active timeslots per 
PRBS code cycle, it will take 8 PRBS cycles (i.e., 8 PRBS 
5 cycles times 4 active t imeslots/PRBS cycle) to induce 32 
Supersteps. 

Similarly, a PRBS code of length 15 has 8 active 
timeslots, and 32 can be evenly divided by 8. For a PRBS 
code of length 31, there are 16 active timeslots, which can 

10 also be evenly divided into 32. 

It should also be noted that, for some of the higher 
flow rates, more than one Superstep of fluid may be delivered 
in each active timeslot. In the example set forth in 
Table A, flow rates over 3.1 ml/hr require two or more 

15 Supersteps per timeslot. The number of Supersteps per 

timeslot is preferably a number that can be evenly divided 
into the total number of Supersteps per pump cycle. 

Where two or more Supersteps of fluid are delivered in a 
single timeslot, the Supersteps are preferably delivered at 

20 the beginning of the timeslot, either as a single bolus of 
fluid or using the "trilling" method discussed below with 
respect to FIG. 13. 

FIG. 7 illustrates a process for generating PRBS codes. 
Each PseudoRandom Binary Sequence (PRBS) consists of 1 and 0 

25 code values, with each 1 or 0 associated with a Timeslot. 
Timeslots associated with code values of 1 produce one or 
more Supersteps of flow at their beginning, while Timeslots 
associated with code values of 0 produce no flow. The period 
of code repetition is typically between 40 and 100 seconds, 

30 varying with the flow rate. The coded flow pattern produces 
a cumulative flow equal to that which would have been 
produced by a nominal constant flow over the time period of " 
the code; 

In a preferred embodiment, PRBS codes are produced only 
35 in integer lengths of M, where M = 2* - 1 with k being an 

integer value. For most applications, codelengths used are M 
=7, 15, or 31. PRBS codes are composed of a series of l's 
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and 0's f with the number of 1's always being one more than 
the number of 0's. The count of l's will always be an even 
multiple of 2. 

In a preferred embodiment, PRBS code sequences are 
5 determined using a feedback algorithm such as that shown in 
FIG. 7. Three key values are used to define the type and 
phase of any PRBS sequence. The codelength M defines the 
length of the PRBS code. The OCTAL generator code determines 
PRBS type, while the Binary Initialization Code (BIC) 
10 determines the phase or starting point within the code. 

Representative codelength M, OCTAL, and BIC values are set 
forth below in Table B. 



TABLE B 



Selected 
Flow Rate 
(ml/hr) 


Codelength 
M 


OCTAL 
Code 


Binary 
Initialization 
Code 
(BIC) 


0.1 - 0.5 


n/a 


n/a 


n/a 


0.6 - 1.4 


7 


13 


111 


1.5 - 3.0 


15 


31 


1001 


3.1 - 6.1 


15 


23 


1010 


6.2 - 12.4 


15 


23 


1101 


12.5 - 24.9 


15 


23 


1100 


25.0 - 49.9 


31 


75 


11101 



The OCTAL values in Table B are from Zeimer and 
25 Peterson, Digital Communications And Spread Spectrum Systems, 
page 390, MacMillan Publishing Co., New York. The OCTAL 
code, which is expressed in. the table in a base 8 (i.e., 
octal) format, determines which of several possible PRBS 
codes of a given length will be formed. 
30 Certain codes and relative rotations (i.e., starting 

points) have been determined to provide optimal performance 
for particular pumps to minimize effects of mechanism and 
sensor generated artifacts. The BIC values in Table B, which 
determine the phase or starting point within the code, were 
35 empirically determined based on flow characteristics of a 
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particular pump, as discussed in greater detail below with 
respect to FIG. 12. 

In the PRBS coding process shown in FIG. 7, the selected 
OCTAL code is converted into a binary format to create an 
5 array of 1 and 0 multipliers. The least significant bit 

(lsb) of the OCTAL code is discarded, and the remaining code 
is arranged, from the next-to-least significant bit (lsb+1) 
(the least significant bit having been discarded) , to the 
most significant bit (msb) in an the Octal Code (binary form) 

10 array 114. The Octal Code (binary form) array is used to 

operate on the outputs of "degree" number of shift registers. 

The Binary Initialization Code (BIC) is used to 
establish the initial value of the shift registers and thus 
determine the first few bits of the PRBS code and the 

15 relationship of the PRBS code to the physical rotation of the 
pump mechanism. The BIC phase value corresponds to the 
motor/rotation monitor being aligned to start Superstep 0. 
As shown in FIG. 7, the BIC is arranged, from least 
significant bit (lsb) to" most significant bit (msb) , in an 

20 array 115. Upon system activation, the BIC array is used to 
preload the generator array elements 116. 

In a first operation, the Octal Code (binary form) array 
elements (114 a to 114 d ) are multiplied with generator array 
elements (116 a to 116 d ) , as shown by multipliers (117 a to 

25 117 d ) . The resulting values are processed through a series 
of EXCLUSIVE/OR ("XOR") function gates 118, with each XOR 
function gate receiving a pair of values as an input and 
outputting either a 1 or a 0 . The XOR gates output a 0 in 
response to an input of a (1 and 1) or a (0 and 0) , and 

30 output a 1 in response to an input of a (1 and 0) or a (0 and 
1) . As shown in FIG. 7 , the two values produced by two of 
the multipliers (117 a and 117 b ) are inputted to one XOR gate 
(118 a ) , and the two values produced by the remaining two 
multipliers (117 c and 117 d ) are initially passed to another 

35 XOR gate (118 b ) . Each o*f these XOR gates (i.e., 118 a and 

118 b ) output a single 1 or 0 value. The outputs of these XOR 
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gates (118 a and 118 b ) are inputted to a third XOR gate 118 c , 
which outputs a single 1 or 0 value. 

The left-most value of the generator array (i.e., the 
value in register 116 a ) is outputted to motor control as the 
5 first element in the PRBS array. The values in all the other 
registers shift one position to the left, leaving the last 
register 116 a empty. The output value from the third XOR 
gate 118 c is then entered into the last register 116 a . 

The process depicted in FIG. 7 is specific to a 4- 
10 element BIC. However, similar processes can be used to 
generate PRBSs using BICs of different lengths. 

The above-cited process is repeated until M elements of 
the PRBS code have been generated. 

The above discussion concentrates on the use of PRBS 
15 codes. However, pseudorandom codes other than PseudoRandom 
Binary Codes are also applicable to the current invention. 
For example, a PseudoRandom code consisting of non-binary 
values, such as a 4 -bit code having values of 2, 4, 0, and 3, 
could be used. Even non-integer values could be used, such 
20 as a 3 -bit code having values of 0.6, 1.2, 0.8. In such an 
embodiment, three boluses of fluid would be delivered, each 
having volumes of proportional to the particular bit (i.e., 
0.6, 1.2, 0.8). Alternatively, the boluses might be of 
approximately the same size, but the spacing between boluses 

2 5 may be varied in proportion (either direct or inverse) with 

the particular bit value . 

A pseudorandom code is one for which its frequency 
response contains energy at all frequencies. This allows for 
a suitable inverse code to be created. Non-binary and non- 
30 integer values are applicable to the current invention, in 
that they can be used to accurately estimate equilibrium 
pressure and determine "virtual" waveforms. However, the use 
of non-binary and non- integer values complicates the 
processing involved. Restricting the codes to PseudoRandom 

3 5 Binary Codes, such as 1101, simplifies processing and 

improves system efficiency. Accordingly, PRBS codes are the 
preferred embodiment of the current invention. 



WO 00/13726 PCMJS99/18334 

-28- " — 

As was shown in FIGS. 3 and 4, flow delivered in 
accordance with a PRBS code results in a corresponding 
pressure waveform, which can be monitored with a transducer 
or other sensor to produces corresponding pressure signals. 
5 The pressure signals can be decoded via a PRBS decoder to 
create a corresponding "virtual" waveform. However, because 
pressure sensors typically provide numerous pressure signals 
per second, hundreds of pressure signals are created for most 
timeslots. Decoding all the pressure signals for each 

10 timeslot is extremely processor- intensive . However, by 
averaging select pressure signals from each timeslot, and 
only decoding a single average value from each timeslot, the 
amount of decoding required can be drastically reduced while 
still maintaining accuracy of the P c . Rather than decoding 

15 each and every pressure signal in a timeslot, an average can 
be taken of just a portion of the timeslot, and this average 
value can be decoded to determine an estimated P G . The 
estimated P D determined using pressure signal averaging is 
remarkably close to the estimate P Q determined by decoding 

20 each and every pressure signal. However, the processing 
requirements are drastically reduced using pressure signal 
averaging . 

FIG. 8 shows a pressure waveform 120 in a preferred 
embodiment of pressure signal averaging, wherein the averaged 
25 portion 122 of the timeslot comprises the final 25% of the 
"tail" portion. This average is known as the Baseline 
Sample. This process minimizes the effects of high frequency 
artifacts, improves effective resolution, and reduces decoder 
operations. 

3 0 It should be noted that the pressure signal averaging 

need not be confined only to 25%, or to just the tail portion 
of the timeslot. Other portions and other percentages of the 
pressure signal could also be averaged and still give 
relatively accurate P c estimates. An accurate P Q estimate 

3 5 can even be made using the average of pressure signals over 
the entire timeslot. When the average is then decoded, the 
resulting estimated P D is remarkably close to the estimated 
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P D determined by averaging only the tail portion of the 
timeslot or by directly decoding each and every pressure 
signal over the entire timeslot. 

For decoding, each averaged value is entered into a 
5 shift register array of length equal to that of the PRBS code 
(i.e., length = M) . The shift register array comprises an 
averaged pressure value from the most recent M timeslots. 
The averaged pressure values are then decoded, and the 
decoded pressure values used to determine an estimated P G . 

10 As an additional enhancement to accuracy of the 

resistance determination, the digital pressure signals 
received from the A/D converter may be premultiplied prior to 
the averaging process. Many A/D convertors have a resolution 
of about only about 0.5 mm Hg per count, which is 

15 insufficient for accurate measurement of low resistances at 
low rates. In a preferred embodiment, additional resolution 
is obtained by premultiplying raw A/D count values by 256 
prior to the averaging process. The premultiplied values are 
then averaged, and the averaged result, now in scaled units, 

2 0 is input to the PRBS decoder to produce scaled-up response 

values. 

FIG. 9 illustrates a process for decoding pressure 
signals using a PRBS decoder. Decoding the resulting 
pressure signals requires M 2 addition or subtraction 

25 operations per resistance calculation, and M array rotate 
operations per resistance calculation. On the first ad 
sample of each timeslot, motor control passes to the decoder 
the current encoding value (0,1). The decoder will enter 
into an M element shift array the value -1 if the encoding 

30 value is 0, or a +1 if the encoding value is 1. The 

resulting array of +l's and -l's, known as the Cbit array, 
will be used for decoding. The Cbit array is initialized to 
all 0's whenever RUN is activated. Since P e estimate 
selection logic will preclude use of the decoded outputs 

3 5 until M timeslots elapse following RUN activation, the Cbit 

array will be filled with correct values by the time it needs 
to be used. 
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The decoded output for an M length array of pressure 
samples P may be expressed mathematically as: 



M - 1 

DecoderOutii) = Y Cbit( f , + M ) P{k) 

' |J modulo M 

k = 0 



where the index i runs from 0 to M-l, and the Oth index 
5 refers to the earliest decoder value. Since the Cbit array 
contains only 1 and -1 values, the multiplication above may 
be replaced by a logical test and an add or subtract 
operation, depending on the value of the Cbit array element. 
FIG. 9 illustrates the decoding process for a simple M=3 

10 code used to decode a 3 -bit array of pressure baseline sample 
averages (BLSAs) . Each sample box contains the index of the 
value to be used in processing, with 0 representing the 
current sample (or unrotated code bit) . The Cbit array 
values are +1, +1, and -1; corresponding to a PRBS code of 1 

15 1 0. 

FIG. 9a shows the alignment of the Cbit array and the 
BLSA array elements within shift registers Cbit and BLSA. 
Each time a timeslot is completed, the Encoder replaces the 
rightmost (current - indexed 0) element of Cbit with its 

20 current value (1 or -1) . Previously held values within Cbit 
are shifted to the left, and the oldest (leftmost) is 
discarded. Similarly, with the completion of each timeslot, 
the Baseline Sample Averager produces a BLSA value which is 
placed in the rightmost (current - indexed 0) element of 

25 BLSA. Previously held values within the BLSA array are 

shifted to the left, and the oldest (i.e., leftmost) value is 
discarded. 

Following the loading and shifting of the arrays, the 
actual decoding of the BLSA array occurs. As shown in FIG. 
30 9a, the product of each corresponding element of the two 
shift register arrays is formed. Then the sum of these 
products is taken. This result will produce the earliest 
"virtual pressure impulse" response, which is indexed as 
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element 0 in the time line graph to the right of the 
registers and denoted Decode (0). 

Referring now to FIG. 9b, the elements of Cbit are 
rotated, so that the newest element is moved into the 
5 position previously occupied by the eldest, and all other 
elements {including the eldest element) are shifted one 
register to the right. All these operations are preferably 
performed with the assistance of a temporary holding register 
(not shown) or equivalent functionality of a microprocessor 

10 so that no data elements are lost. 

Again, the corresponding products for each shift 
register position are formed, and the sum of the products 
computed. This produces the next-to-earliest "virtual 
pressure impulse" response element, which is indexed as 

15 element 1 in the time line graph to the right of the 
registers and denoted Decoded). 

The process of rotation of the elements, multiplication 
of corresponding elements, and summing the products is 
repeated for the remaining Cbit elements, which in the 3- 

20 element array of FIG. 9 is only one additional element. As 
shown in FIG. 9c, the resulting element is indexed as element 
2 in the time line graph to the right of the registers and 
denoted Decode (2) . 

The above-cited decoding process is known generically as 

25 "circular convolution," and is well known in the signal 
processing art . 

Another example of PRBS decoding is shown in Table C. 
In the example shown, a PRBS code of 1110100, which was 
used to generate a flow in a conduit, is used to decode an 

30 array of measured pressure signals, which in the example 
shown is (4 5 6 4 5 4 3). For the first row, i.e., row A, 
the binary values (1X10 10 0) are replaced with a plus 
(+) or minus (-) sign, with a plus ( + ) sign replacing 1 
values, and a (-) sign replacing 0 values, resulting in the 

35 series (+ + + - + --). " For each subsequent row (i.e., rows 
B through G) , the resulting series of +/- signs are each 



i 
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shifted one place to the right (with the +/- sign in the 

rightmost column transferring to the leftmost column) . 

TABLE C 



10 
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+ 
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+ 
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+ 




+ 








1 



The resulting +/- signs in rows A through G are then 
applied to the pressure values to be decoded, and the 
resulting values are added together across each row (A 
through G) . Thus, row A, when applied to an array of seven 

15 pressure averages pressure values (4, 5, 6, 4, 5 , 4, 3), 

translates to (+4, +5, +6, -4, +5, -4, -3), which when added 
together results in a value of +9. Row B translates to (-4, 
+5, +6, +4, -5, +4, -3), which when added together results in 
a value of +7. Row C translates to (-4, -5, +6, +4, +5, -4, 

20 +3), which when added together results in a value of +5. 
Rows D through G result in values of 3, 3, 3, and 1, 
respectively. Accordingly, the decoded array, which 
represents the decoded pressure waveform, is (9753331). 
(Note that the average pressure values actually shown in 

25 Table C are not based upon actual measurements, but are 
sample numbers that are used merely for illustrative 
purposes . ) 

After the pressure signals have been decoded, they can 
be used to determine an estimated P G . Various methods can be 
30 used to estimate P 0 , with the preferred method varying 
depending on the user-selected flow rate, the time since 
start up, the values of the decoded pressure responses, and 
previous resistance estimates. In a preferred embodiment of 
the invention, the method for determining estimated P is 
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selected as a function of various parameters, as set forth in 
Table D. 

Various methods for determining the estimated P D are 
discussed below. The logic for selecting a particular method 
5 for determining P Q is shown in Table D. 

For user selected flow ranges of 0 . 5 ml/hr or less, 
there is typically sufficient time in each timeslot for the 
monitored pressure to decay to P D . However, transient 
pressure and sensor anomalies can cause brief spikes or other 

10 inconsistencies in the monitored pressure. To reduce the 
impact of such monitored pressure inconsistencies, the 
estimated P 0 is determined as the average of the final 
portion, such as the final 25%, of the pressure readings in 
the timeslot. No pseudorandom coding or decoding is 

15 involved. Determining P G for low flow rates, such as flow 
rates of 0.5 ml/hr or less, is discussed* in greater detail 
below with respect to FIGS. 14 and 15. 

Another group of techniques for determining P Q that is 
useful when pseudorandom coding and decoding are used is 

20 referred to as Three Point Block Averaging. Three point 
block averaging is a process whereby the pressure samples 
that are output from the decoder, which has M number of 
samples (i.e., a length of M) , is divided into three blocks 
of samples, and the average of each block of samples is 

25 taken. The resulting block averages are then used to 
determine an estimated P c . 

Of course, M is not always evenly divisible by three. 
However, where M is not evenly divisible by three, M-l will 
always be divisible by three. Accordingly, where M cannot be 

30 evenly divided by three, the last sample of the decoded 
output is discarded, which gives a remaining number of 
samples that is evenly divisible into three blocks. 

As shown in FIG. 10, the average of each block of an 
output of decoded samples 128 is shown, where M = 15. The 

35 resulting averages are defined as follows: 
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A 0 = the average of the first block of samples; 
A x = the average of the second block of samples; 
A 2 = the average of the third block of samples. 
Intermediate values (i.e., the differences between adjacent 
> average values) are as follows: 

t>i = A 0 - h x 
D 2 = A a - A 2 

Using all or portions of the above values, i.e., A 0 , A lt 
A 2 , D 1# and D 2 , P Q can be estimated. Various techniques for 
using A 0 , A w A 2 , D 1# and D 2 to determine the estimated P Q are 
set forth below. Which technique to use depends on the 
particular parameters, such as pump characteristics, flow 
rate, NSAD status, etc. In a preferred embodiment of the 
invention, the technique for determining P G is selected using 
15 the parameters set forth in Table D. 

Formulas for determining the estimated P D that use three 
point block averaging include the following: 



10 



20 



(A>) 



2 



D 2 - D x 



P Q = Average of A, and A 



Another method for determining P G , but without using 
three point block averaging, involves averaging just the tail 
portion of the pressure signals from each timeslot. These 
tail averages are not decoded. Instead, the tail averages 
are themselves averaged, resulting in a single average value. 
The average value is then used as the estimated P c . This 
method is generally effective under conditions where low 
25 resistances are expected, such as where the most recent 
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previous measurement from the output of the resistance median 
filter indicated a low resistance value. 

As shown by the above, there are a variety of methods to 
estimate P 0 . The suitability of particular methods depends 
5 on various conditions, such as flow rate, set compliance, 
etc. Determining the optimum method can require extensive 
testing of the particular equipment. Table D shows various 
preferred methods for estimating P D , with the preferred 
method a function of selected flow rate, NSAD factor (i.e., 

10 stability of the decoded signal) , calculated resistance, and 
time the system has been in operation. 

When large changes in baseline or impedance occur (due 
to such factors as noise, large offset changes, or hard 
occlusions), the decoded output can become unstable. To 

15 prevent gross miscalculations of resistance under such 

conditions, an alternate estimation method must be selected. 
A preferred embodiment of a P Q estimation selection logic, 
such as the one shown in Table D, requires a method to 
determine stability of the decoded signal. This is achieved 

20 by computing an "estimate noise" parameter known as 

Normalized Sum Absolute Delta's (NSAD) . NSAD is defined as 
the sum of the absolute value of the differences D x and D 2 
divided by the product of codelength (M) and the number of 
Supersteps per timeslot (to normalize for different codes and 

25 drive patterns) . NSAD is thus defined by the following 
formula : 

NSAD = ( I J, I * \D 2 \ ) nrmHg 

CodeLength(M) * Super stepsPer Time 



NSAD is compared to a threshold NSAD threah/ which is itself 
dependent on the selected flow rate. Preferred embodiments 
of values of NSAD threflh for various selected flow rates are 
30 shown in Table E. If NSAD exceeds NSAD thre8h , a "true" 

condition is input into the P^ determination method selection 
logic process, as shown in Table D. 
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TABLE E 



Selected Flow Rate 
(ml/hr) 




0.1 - 0.5 


n/a 


0.6 - 1.4 


4 


1.5 - 3.0 


4 


3.1 - 6.1 


7 


6.2 - 12.4 


7 


12.5 - 24.9 


7 


25.0 - 49.9 


7 



When fluid resistances are very low, the difference 
values Dj and D 2 become very small. Where P D is determined 
using the formula with the value (D 2 - DJ in the denominator 
This can cause the difference D 2 - D x to approach zero. 

15 Where that difference (i.e., D 2 - D x ) is used in the 
denominator of an equation for determining P D , such as 
Equation 8 above, a large and unstable P G estimate can 
result. To identify situations that might yield such 
unstable results, the NSAD value is compared with a minimum 

20 NSAD threshold value NSAD rainthresh , which in a preferred 

embodiment is equal to 0.25 mm Hg/ (M * SuperstepPerTimeslot ) , 
and the logical result is used in the P c estimation method 
selection processor outlined in Table D. 

When the resistance becomes relatively low, the most 

25 stable method for estimating P G uses a direct average of the 
M available Baseline Sample Average (BLSA) values as the 
estimated P e , without any decoding of the pressure averages. 
This technique is used when output from the Resistance 
Estimation Filter is less than a resistance threshold Rt hreS h/ 

30 which in a preferred embodiment is 600 fluid ohms. 

Following a period of no flow, such as just after system 
start up, the output from the PRBS decoder will be incomplete 
and inaccurate until at least M timeslots have elapsed since 
the mechanism was started. Accordingly, there is 

35 insufficient data in the system to accurately calculate an 

estimated P c . Where less than M timeslots have elapsed since 
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the mechanism was started, other methods for selecting P D 
estimation value must be used. Two such methods are set 
forth in Table D. In one embodiment, for selected flow rates 
between 0.6 to 50 ml/hr, the P e value used is the pressure 
5 that was measured in the conduit just prior to the beginning 
of pump operation (i.e., just before system run). For rates 
of 0.5 ml/hr or less, the Baseline Sample Average for the 
individual timeslot alone is used, in combination with the 
special "low flow rate" mode of resistance calculation 
10 discussed in greater detail below with respect to FIGS. 14 
and 15. 

Spurious signals and other transient anomalies can 
occasionally result in an inaccurate estimated P Q for 
occasional timeslots. As shown in FIG. 11, the impact of 

15 such inaccurate P D estimates can be reduced and sometimes 
eliminated by using a median filter. The median filter 
includes an array for storing the most recent P Q values. In 
a preferred embodiment shown in Table F, the P D filter array 
has a size M, so that it stores the P G values for one 

2 0 codelength. 



TABLE F 



Selected Flow Rate 
(ml/hr) 


Filter Array 
Length 


0.1 - 0.5 


n/a 


0.6 - 1.4 


M (7) 


1.5 - 3.0 


M (15) 


3.1 - 6.1 


M (15) 


6.2 - 12.4 


' M (15) 


12.5 - 24.9 


M (15) 


25.0 - 49.9 


M (31) 



The median filter process shown in FIG. 11 involves an 
array of length 7. The array stores the most recent 7 
values. The relative age of individual data elements is 
maintained in a linked list, shown at (A) , that shows the 
35 relative ages from oldest to newest. 
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The data elements are maintained and sorted in an array, 
shown at (B) , from smallest to largest. When a new value is 
presented, such as the value "11" shown at (C) , the oldest 
value is discarded, which in the example shown at (D) is the 
5 value "25". The new value takes the place of the discarded 
old value. Note that in the example shown, the array was 
preloaded before operation with the value "25" in each array 
slot. Prior to introduction of the most recent value (i.e., 
"11"), only four new values ("9", "12", "18", and "43") were 

10 added to the array. 

The array with the newly added value ("11") is sorted, 
such as through a typical bubble sorting method, to place all 
values in order from smallest to largest, as shown at (E) . 
The resulting ordered array, shown at (F) , has all elements 

15 from smallest to largest. The center value, which in the 
example shown is "18", is the new median filtered value, as 
shown at (G) . 

Note that filtering could also be performed by other 
methods, such as simple averaging of the current P c value 

20 with other recent P Q values. However, median filtering has 
certain advantages over averaging, particularly in 
eliminating the effects of extreme values, such as the "43" 
that appeared in the example shown in FIG. 11. 

Another method for improving the accuracy of the P 0 and 

25 resistance calculations involves the use of the BIC code, 

which is used to reduce the impact of revolution synchronous 
noise sources . Revolution synchronous noise sources are 
anomalies in pressure signals that occur predictably within 
specific Supersteps. These noise sources can corrupt the 

3 0 decoded waveform, causing large inaccuracies in estimated P c . 
Revolution synchronous noise sources include mechanism- 
induced sensor offset shifts and variations in flow between 
Supersteps . 

Variations in flow volume Qss between Supersteps can 
35 range from small variations to situations where a Supers tep 
actually produces negative flow. Although the Supersteps are 
selected to deliver approximately equal volumes of fluid, as 
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shown in FIG. 12a, fluid delivery varies between Supersteps 
over time, depending on various system parameters such as 
downstream pressure. In a preferred embodiment, the 
Supersteps are predetermined based upon pump characteristics, 
5 and are not dynamically adjusted to compensate for time- 
variant changes in fluid delivery volumes. However, under 
some circumstances, such as where the downstream pressure is 
large, the volumes delivered by different Supersteps can vary 
widely, even to the point where one or more Supersteps may 

10 deliver negative volumes of fluid. For example, while a pump 
operating with relatively low downstream pressures may 
produce the substantially equal fluid delivery volumes shown 
in FIG. 12a, that same pump operating with higher downstream 
pressures may deliver varying fluid delivery volumes such as 

15 those shown in FIG. 12b. 

Although most Supersteps either experience no change or 
only a very small change in delivered volumes due to 
downstream pressure, one or rrtore Supersteps can see 
variations in their delivered volume. For example, in the 

20 hypothetical 200-step, 32-Superstep pump cited by way of 

example in FIGS. 12a and 12b, Superstep 29 (SS 29 ) can go from 
a small positive delivered volume Qss 29 at low downstream 
pressures (FIG, 12a) to a large negative delivered volume 
Qss 29 at higher downstream pressures (FIG. 12b) . Such time 

25 and pressure-variant changes, i.e., "artifacts, " in Superstep 
delivered volumes can cause problems with accurately 
determining equilibrium pressure P Q and resistance R. 

Noise synchronous sources are generally predictable in 
nature. The problems caused by revolution synchronous noise 

30 sources with determining P 0 and R can largely be eliminated 
by pairing a "noisy" Superstep (such as the negative-volume 
Superstep SS 39 from FIG. 12b) with a particular PRBS 
codestep: This is accomplished by the Binary Initialization 
Code ( M BIC n ), which is used to establish the initial value of 

35 the shift registers and thus determine the relationship of 

the code to the physical rotation of the pump mechanism. The 
BIC phase value typically corresponds to the motor/rotation 
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monitor being aligned to start Superstep 0 (SS 0 ) . Correct 
phasing of the code relative to the mechanism helps to 
minimize artifact-induced errors in the P Q estimate. In a 
preferred embodiment, the BIC is empirically determined for 
5 particular pumps and flow rates, depending on the noise 

synchronous sources for the particular system. For the 200- 
step IVAC Signature Edition Pump (versions 7100 and 7200)., 
preferred BICs for various flow rates were shown previously 
in Table B. 

10 In order to maximize the time allowed for the pressure 

to decay while still maintaining timeslots of a reasonably 
small size, it is generally desirable to deliver fluid in a 
narrow time range, such as the very beginning, of a timeslot. 
Such delivery preferably occurs by delivering one or more 

15 Superstep volumes (Qss) using a single 

acceleration/deceleration current waveform. For example, 
FIG. 13a shows four Supersteps delivered as a single 
acceleration/deceleration waveform 130. The use of 
acceleration/deceleration current waveforms is described in 

20 detail in pending U.S. patent application Serial No. 
08/526,468, which is incorporated herein by reference. 

However, where more than one Superstep occurs in a 
single timeslot, or where a single Superstep involves the 
delivery of a substantial amount of fluid, delivering such a 

25 large amount of fluid in a single acceleration/deceleration 
waveform 130a, as shown in FIG. 13a, may create undesirably 
high instantaneous flow rates, which can create "sympathetic 
flow" from the primary fluid source during secondary 
delivery. Additionally, using a single 

30 acceleration/deceleration waveform can create undesirably 
high instantaneous pump speeds. 

To reduce instantaneous pump and fluid injection speeds, 
the fluid delivery may be "trilled." Trilling involves 
breaking the single waveform into a plurality of smaller 

35 waveforms 130b # as shown in FIG. 13b, with each small 

waveform preferably following an acceleration/deceleration 
waveform. In the embodiment shown in FIG. 13b, four 
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Superstep volumes (4 Qss) , which were delivered as a single 
bolus waveform 130a in FIG. 13a, are delivered as four 
separate waveforms 130b, with each small waveform 130b 
delivering one Superstep volume (Qss) of fluid. The 
5 individual waveforms 130b are preferably spaced immediately 
adjacent to each other so that fluid delivery occurs in a 
narrow time period, with the remainder of the timeslot 
allowing for the pressure to settle/decay to allow for more 
accurate estimation of P 0 . 

10 In addition to an accurate P Q value, determining 

resistance also requires determination of the sum of the 
pressure signals (i.e., pressure samples) for each timeslot. 
Summing the pressure signals does not require the signals to 
be decoded accordingly, the undecoded, or "raw," signals 

15 are used for the summation process. During each timeslot, 
the sum of all or a substantial portion of the pressure 
samples collected during that timeslot is computed. The sum 
value for the timeslot is entered as one element in a length 
M array, known as the TimeSlotSumArray (TSSA) . In a - 

20 preferred embodiment, this array is initially set to all 0's 
to ensure well behaved startup. For each resistance 
computation, the sum of TSSA (i.e., the sum of all pressure 
samples over the past M timeslots) is updated and used for 
the resistance computation. 

25 FIG. 14 shows a sample TSSA having a length of 7, 

corresponding to a PRBS codelength (M) of 7, used to 
determine a pressure summation for a corresponding undecoded 
pressure waveform 132. Each timeslot has a plurality of 
undecoded pressure samples 134 that define the pressure 

30 waveform. Each element 136 of the TSSA is used to store the 
sum of the undecoded pressure samples 134 from a particular 
timeslot. All elements of the TSSA are then added together 
to determine an pressure summation 138. 

Note that the TSSA is preferably used to keep a moving 

35 or running total of the pressure sums for the most recent 
timeslots. In other words, the TSSA is updated each time 
that pressure signals from a new timeslot are received, with 
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the data from the "oldest" timeslot being replaced by the new 
data . 

Once the pressure summation and P c are known, resistance 
can be computed as: 

5 

277.77 * Y^ TimeSlotSumAlx:a yj ~ P Q est*Yl SampleSuirtj 

SamplesPerSec * Supers tepsPer Timeslot * 

Super step 



where: R = resistance expressed in fluid ohms, 

P = pressure in mm Hg, and 

277.77 converts mm Hg sec//xliter to fluid ohms 
(i.e., mm Hg hour/liter). 

10 It is critical that the value of ESampleSum^ be 

identical to the total number of pressure samples 
contributing to ETimeSlotSumArrayj 

For each timeslot, a new resistance estimate is 
calculated. However, spurious signals and other transient 

15 anomalies between timeslots can cause inaccuracies in 

calculated resistance values. To reduce the impact of these 
anomalies, the calculated resistance value can be processed 
with a median filter. Accordingly, each new resistance 
estimate is entered into a resistance median filter array. 

20 The resistance median filter operates similarly to the P 0 

median filter, producing a filtered resistance value that is 
generally an accurate estimate of fluid resistance in the 
system. As was the case with the P D filter, the median 
filter determines the median of the resistance estimates in 

25 the array, and this median becomes the filtered resistance 
value, the resistance median filter array has a length that 
is rate dependent, with ^values for a preferred embodiment set 
forth in Table G. The array elements of the filter are 
preferably initially set to 0 to minimize startup transients. 

30 When the user selected rate is changed without stopping the 
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fluid injection system, all filter array elements are 
initially set with the most recent filtered resistance value. 
In a preferred embodiment, the filtered resistance value is 
presented on the display and alarm processing systems. 



5 TABLE G 



Selected Flow Rate 
(ml/hr) 


Resistance 
Filter Array 
Length 


0.1 - 0.5 


4M (16) 


0.6 - 1.4 


4M (28) 


1.5 - 3.0 


4M (60) 


3.1 - 6.1 


3M (45) 


6.2 - 12.4 


3M (45) 


12.5 - 24.9 


3M (45) 


25.0 - 49.9 


3M (93) 



15 In addition to determining resistance, the system can 

also determine other flow parameters, including compliance, 
system impedance, etc. For example, the system could be used 
to determine a Time Constant TC (which equals Resistance 
times Compliance) using the following formula: 

f t (P(t) - P Q ) dt 
TC = ± - 

j (Pit) - P Q ) dt 



20 The above reduces to the discrete form: 

H 

£ k AT (P(k) - P Q ) 
TC = i^-i 

M 
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where : 

k = the index range; 

AT = the timeslot interval (i.e., timeslot length); 
M = the number of samples; 
5 P(k) = the M samples from the decoder; and 

P Q = estimated equilibrium pressure. 
Using the above determination of TC in combination with the 
separately determined resistance (such as the filtered 
resistance value) , system compliance can be easily 

10 determined. Because TC = Compliance times Resistance, 
Compliance can be simply computed as TC divided by R. 

An additional benefit of using a pseudorandom code is 
the ability to discriminate pressure signals caused by other 
sources of noise, even pressure signals caused by other pumps 

15 operating on a common fluid infusion system. In the 

embodiment depicted in FIG. 17, a fluid infusion system or 
assembly 150 includes two fluid infusion segments 152a, 152b, 
each of which comprises a separate fluid source 22a, 22b, 
separate processors 30a, 30b, separate pump motors 14a, 14b, 

2 0 and separate pump mechanisms 16a, 16b acting on separate 
fluid lines 12a and 12b. The separate fluid lines 12a and 
12b feed into a common fluid line 12c for infusion to a 
patient 24 through a common cannula 26. Accordingly, the 
fluid lines of the two systems essentially form a common 

25 fluid conduit. Thus, fluid flow induced by pump mechanism 
16a causes corresponding pressure responses that are 
monitored not only at pressure sensor 34a in conduit 12a but 
also at pressure sensor 34b in conduit 12b. Similarly, fluid 
flow induced by pump mechanism 16b causes corresponding 

30 pressure responses that are monitored not only at pressure 
sensor 34b in conduit 12b but also at pressure sensor 34a in 
conduit 12a. Without proper filtering, such pressure 
response crosstalk can degrade the accuracy of resistance and 
other measurements . 

35 Pseudorandom coding and decoding inherently filters out 

such crosstalk, making it possible for each separate infusion 
segment 152a, 152b to accurately determine flow resistance. 
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For example, infusion segment 152a preferably operates using 
a pseudorandom code that causes a resulting fluid flow 
pattern that is distinct from the flow pattern induced by 
infusion segment 152b. By decoding the resulting pressure 
5 response in accordance with the pseudorandom code, the 

processor 3 0a of infusion segment 152a inherently filters out 
much of the crosstalk caused by infusion segment 152b. 
Accordingly, resistance can be determined with relative 
accuracy. Accuracy can be further enhanced by applying the 

10 median filter technique shown in FIG. 11. Thus, the 
individual processor 30a, 30b of each of the infusion 
segments 152a, 152b can individually determine resistance by 
filtering the monitored pressure signals with the 
pseudorandom code used by the particular processor 30a, 30b 

15 to generate fluid flow with the respective pump mechanism 
16a, 16b. 

FIG. 17 shows both of the infusion segments 152a, 152b 
as having processor-controlled infusion mechanisms. In one 
embodiment , each processor 30a, 30b controls fluid flow in 

20 accordance with a pseudorandom code. This further enhances 
the ability of each processor to filter out crosstalk from 
the opposing infusion segment. However, if both processors 
30a, 30b were using the same pseudorandom code, the ability 
of the coding/decoding process to filter out crosstalk would 

25 be seriously degraded. Additionally, there are certain 

pairings of pseudorandom codes that, when each is used in a 
separate infusion segment (such as 152a, 152b) , can degrade 
or enhance the ability of the coding/decoding process to 
filter out crosstalk. Accordingly, in a further embodiment, 

30 the processors 30a, 30b cooperate to ensure that each 

infusion segment 152a, 152b is infusing fluid in accordance 
with a different pseudorandom code than the other infusion 
segment. The processors may further cooperate to select 
pseudorandom codes that enhance each processors ability to 

35 filter out the crosstalk created by the other infusion 
segment. For example, the processors 30a, 30b may 
purposefully select different Octal codes that generate PRBS 
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codes which enhance the ability of the processors 30a, 30b to 
filter out the crosstalk via the coding/decoding process. In 
another embodiment, the separate infusion segments 152a, 152b 
may be controlled by a single processor the coordinates the 
5 pseudorandom coding used to drive each pump mechanism 16a, 
16b. 

In the embodiment depicted in FIG. 17, both infusions 
segments 152a, 152b include a processor 30a, 30b that 
controls fluid infusion. However, resistance can be 

10 determined where only one of the infusion segments includes a 
processor that controls fluid infusion according to a 
pseudorandom code. For example, if infusion segment 152a was 
a conventional infusion device that did not operate according 
to a pseudorandom code, infusion segment 152b could still 

15 determine a resistance so long as processor 30b used 
pseudorandom coding and decoding. 

The pressure sensor used by a particular processor (such 
as 30b) does not necessarily have to be positioned in the 
immediately adjacent fluid conduit (such as conduit 12b) . 

20 Because pressure responses are often similar throughout the 
fluid infusion system, the pressure sensor might be 
positioned almost anywhere in the infusion system that is 
downstream of the pump mechanisms 16a, 16b. For example, the 
pressure sensor might be positioned in common conduit 12c, as 

25 shown by sensor 34c. Thus, processor 3 0b could receive the 
pressure signals provided by sensor 34c, sum and decode those 
signals, and accurately determine resistance and other 
parameters of the system. 

While FIG. 17 shows two fluid infusion segments 152a, 

30 152b, pseudorandom coding (such as PRBS coding) could be used 
in infusion systems having. three or more fluid infusion 
segments. Although each additional infusion segment adds 
additional "noise" or crosstalk to the system, the 
pseudorandom coding/decoding process can eliminate much of 

35 the crosstalk. 

The above-cited system and method for medium flow rates, 
whereby pseudorandom coding and decoding is used with the 
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monitored pressure response, is effective for determining 
fluid resistances and other parameters over a wide range of 
flow rates. However, for very low flow rates, the pressure 
produced is quite small. Consequently, small errors in the 
. 5 P Q estimation can produce large errors in the computed 

resistance. At extremely low rates, such as those below 0.5 
ml/hr, mechanical -coupling induced errors and thermal drift 
of the pressure transducer can produce very significant 
errors in a PRBS encoded output, leading to large errors in 
10 the computed resistance. Accordingly, an alternate means for 
resistance computation may be preferable. 

FIG. 15 illustrates a system and method in a preferred 
embodiment of the invention for determining resistance for 
very low flow rates, wherein resistance determination is 
15 performed on a timeslot by timeslot basis. The drive 

controller does not use PRBS encoding. Instead, timeslot 
length is maximized, all timeslot lengths are set to equal 
values, and each timeslot contains a single Superstep. This 
allows continually updated resistance estimates to be 
20 produced with each timeslot. 

In the system for infusing fluid to a patient 24 shown 
in FIG. 15, the user selects a flow rate at a keyboard or 
other rate input device 28. The selected rate is then used 
by the motor controller 46 to determine timeslot length, 
25 which is preferably maximized. The motor pump 14 causes the 
pump mechanism 16 to act on the conduit 12 to deliver one 
Superstep of fluid in each timeslot, preferably as a single 
bolus toward the beginning of the timeslot. By maximizing 
timeslot length, sufficient time is permitted for the 
30 monitored pressure to decay to the equilibrium pressure. 

The sensor 34 monitors the resulting pressure, and 
provides pressure signals that are A/D converted 36 and used 
for pressure summation 48 and tail averaging 54. Tail 
averaging is performed only over a portion of the timeslot, 
35 as shown in FIG. 16. 

FIG. 16 illustrates a sample pressure waveform 140 
created using the system and method shown in FIG. 15. The 



WO 00/13726 PCT/US99/18334 

-49- — 

resulting pressure waveform 14 0 is shown for three subsequent 
timeslots TS n . lf TS n , and TS^ . Each timeslot shows an 
initially increasing pressure wave 142 that decays down to an 
equilibrium pressure. The last portion of the decaying 
5 pressure waveform in the timeslot is known as the tail 144, 
so that an average taken over this last portion is also known 
as the Timeslot Tail Average (TTA) . The equilibrium pressure 
P Q is determined as an average of pressures in order to 
prevent transient pressure anomalies, such as those caused by 

10 patient movement, from skewing the determined equilibrium 
pressure. In the embodiment shown in FIG. 16, the TTA for 
each timeslot is determined as the average of the monitored 
pressures over the last 25% of the timeslot. 

Referring again to FIG. 15, the TTA is calculated 54 and 

15 used in combination with the full sum of pressures over the 
entire timeslot in fluid resistance calculation 50. In a 
preferred embodiment, the resistance is determined using the 
following formula: 

Sampl es/ Timesl o t-1 

277.77 * 52 Pressure] -TTA * Samples/ Timeslot 

R = ill 

SamplesPerSec * ^ 

Superstep 



20 where: R = resistance expressed in fluid ohms, 

P = pressure in mm Hg, and 

277.77 converts mm Hg sec//xliter to fluid ohms 
(i.e., mm Hg hour/liter) . 
This technique for determining fluid resistance is generally 
25 accurate in detecting soft occlusions and other abnormalities 
in the fluid delivery system. However, in the case of a hard 
downstream occlusion, the TTA as well as the other pressure 
sample values of a timeslot would be elevated by the full 
amount of pressure developed by the pump's attempt to 
30 displace the Superstep of fluid volume into the occluded 
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fluid delivery system. Accordingly, the pressure waveform is 
often substantially the same between adjacent timeslots, 
except that the entire waveform is elevated. This is shown 
in FIG. 16, where the pressure waveform portion in the first 
5 timeslot (i.e., TS^) is substantially the same as the 

pressure waveforms in the second and third timeslots (i.e., 
TS n and TS n+i ) except for the waveform in each timeslot being 
slightly elevated over the previous timeslot. Thus, use of 
just the TTA alone might prevent detection of the occlusion. 

10 To remedy this problem, a threshold value is established 

to detect relatively large changes in pressure values between 
adjacent timeslots. In a preferred embodiment, a threshold 
value or slew limit is determined by adding a selected value 
to the TTA (i.e., the P c ) from the previous timeslot. If the 

15 TTA from a timeslot exceeds the slew limit threshold (i.e., 

* 

exceeds the TTA of the previous timeslot by more than the 
selected value), a hard occlusion may be indicated. This 
technique allows for accurate detection of hard occlusions 
while tracking slow drift (i.e., pressure changes of less 
20 than 1 mm Hg/hr which occur due to very slow processes, such 
as gradual inclusion of moisture, curing of adhesive bonds, 
migration of implanted ions within the silicon structure, 
etc.) The technique also avoids errors due to mechanism 
coupling . 

25 In the embodiment shown in FIG. 16, the selected value 

added to the previous Timeslot Tail Average to set the slew 
limit threshold is 2 mm Hg. Pressure increase between 
adjacent timeslots due to a hard occlusion will generally be 
much larger than 2 mm Hg, while the magnitude of Superstep to 

3 0 Superstep pressure errors due to mechanical coupling and 

similar factors is typically less than 2 mm Hg. Accordingly, 
the system compares the TTA of a timeslot to the TTA of the 
previous timeslot. If there is an increase of 2 mm Hg or 
more in the TTA of the more recent timeslot relative to the 

35 immediately previous timeslot, then the current P Q value will 
be limited to the previous TTA + 2 mm Hg. Since hard 
occlusions produce significantly larger pressure increases, 
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the resulting resistance computed from this limited P Q 
estimate will be very high, indicating a hard occlusion. 

By way of example, TTA„ of timeslot TS n in FIG. 16 is 
only 0.5 mm Hg greater than TTA^i of timeslot TS n . t . Such an 
5 increase is likely due to various error factors, and is not 
necessarily indicative of a hard occlusion. However, TTA^ 
of timeslot TS n+l is 2.5 mm Hg higher than TTA,, of timeslot 
TS n . Because TTA^ is more than 2 mm Hg greater than TTA„, 
the slew limit has been exceeded. As shown in FIG. 15, if 

10 slew limit examination 146 determines that there is too great 
a difference between adjacent TTA values, the current P Q 
value will be limited to the previous TTA + 2 mm Hg. When 
the limited P Q value is used to determine resistance, a large 
resistance value will likely result, and the system will 

15 activate an occlusion alarm. 

Although preferred and alternative embodiments of the 
invention have been described and illustrated, the invention 
is susceptible to modifications and adaptations within the 
ability of those skilled in the art and without the exercise 

20 of inventive faculty. Thus, it should be understood that 
various changes in form, detail, and usage of the present 
invention may be made without departing from the spirit and 
scope of the invention. Accordingly, it is not intended that 
the invention be limited, except as by the appended claims. 



10 
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WHAT IS CLAIMED IS: 

1. A system for monitoring at least one flow parameter 
in a fluid delivery assembly in which a flow control device 
acts on a. fluid conduit to control the flow of fluid through 
the conduit, the system comprising: 

a pressure sensor coupled to the conduit, for providing 
pressure signals in response to the pressure sensed in the 
conduit; and 

a processor that controls the flow control device to 
cause flow in accordance with a pseudorandom pattern, 
receives the pressure signals, and processes those pressure 
signals to determine a flow parameter. 



2. A method for determining at least one flow 
parameter in a fluid delivery assembly in which a flow 
control device acts on a fluid conduit to control the flow of 
fluid through the conduit, the method comprising the steps 

5 of: 

(a) controlling the flow control device to cause flow 
in accordance with a pseudorandom pattern; 

(b) monitoring the pressure in the conduit and 
generating signals in response to the monitored 

10 pressure; 

(c) processing the pressure signals to determine a flow 
parameter. 
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